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Abstract 

, We examine the phenomenological imphcations at cohiders for the existence of 
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, higher-derivative gravity terms as extensions to the RandaU-Sundrum model. Such 

OO 

' terms are expected to arise on rather general grounds, e.g., from string theory. In 
CN ' 

' 5-d, if we demand that the theory be unitary and ghost free, these new contributions 

■ to the bulk action are uniquely of the Gauss-Bonnet form. We demonstrate that the 
usual expectations for the production cross section and detailed properties of graviton 

Oh: 

D . Kaluza-Klein resonances and TeV-scale black holes can be substantially altered by 



existence of these additional contributions. It is shown that measurements at future 
^ ■ colliders will be highly sensitive to the presence of such terms. 
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1 Introduction 



The Randall-Sundrum(RS) modelpP provides a geometric solution to the hierarchy problem 
through an exponential warp factor whose magnitude is controlled by the separation, vrrc, 
of two 3-branes embedded in 5-d Anti-deSitter space, AdS^. It has been shown that this 
interbrane distance can be naturally stabilized at a value necessary to produced the exper- 
imentally observed ratio of the weak and Planck scales|2]. In its simplest form, SM matter 
in the RS model is confined to one of the 3-branes while gravity is allowed to propagate 
in the bulk. A very generic signature of this kind of scenario is the existence of TeV-scale 
Kaluza-Klein(KK) excitations of the graviton with inverse TeV-scale couplings to the SM 
fields. These states will appear as a series of spin-2 resonances in a number of processes 
that should be observable at both hadron and e^e^ colliders which probe the TeV-scale. 
The masses and couplings of these KK graviton states will be determined with reasonably 
high precision by future collider measurements. Another possible RS signature is the copious 
production of TeV scale black holes, though this is not a unique feature of the RS modelj^]. 

One can easily imagine that this simple RS scenario is incomplete from either a top- 
down or bottom- up perspective. We would generally expect some 'soft' modifications to 
the details of the picture presented above, hopefully without disturbing the nice qualitative 
features of the model. One such extension of the basic RS model is the existence of higher 
curvature terms which might be expected on general grounds from string theory or other 
possible high-scale completions. Once we open the door to such possibilities the number of 
potential new terms in the action can grow rather rapidly as the number of scalars that 
we can form from products of the curvature and Ricci tensors as well as the Ricci scalar 
are enormous. A certain general class of such invariants with very interesting properties 
was first generally described by Lovelockjn] and, hence, are termed Lovelock invariants. 
(They are also sometimes referred to, apart from a factor of ^/—g, as Euler densities since 
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their volume integrals are related to the Euler characteristics.) The Lovelock invariants 
come in fixed order, n, which we denote as that describes the number of powers of 
the curvature tensor, contracted in various ways, out of which they are constructed. Given 
a space of dimension D the order of these invariants is constrained: For D = 2n, the 
Lovelock invariant is a topological one and leads to a total derivative in the action whereas 
all higher order invariants, D < 2n — 1, can be shown to vanish identically by various 
curvature tensor index symmetry properties. For D > 2n + 1, the £„ are dynamical objects 
that once introduced into the action for gravity can be shown to lead only to second order 
equations of motion as is the case for ordinary Einstein gravity, i.e., no terms with derivatives 
higher than second will appear in the equations of motion due to their presence. Generally, 
arbitrary invariants formed from ever higher powers of the curvature tensor will lead to 
equations of motion of ever higher order, i.e., ever more co-ordinate derivatives of the metric 
tensor and graviton field, e.g., terms with quartic derivatives. Such theories will lead to 
very serious problems with both the presence of ghosts as well as with unitarityj^. The 
Lovelock invariants are constructed in such a way as to be free of these problems making 
them very special and are found to be just the forms taken by the higher order curvature 
terms generated in perturbative string theory|31|S]. This is just what we may have expected 
if string theories are to avoid these unitarity issues. Indeed, if one tracks potential ghost 
generating contributions through the fog of higher derivative terms in the Lovelock invariants 
one sees that they vanish identically. 

In the case of 4-d, apart from numerical factors, Cq = 1 while Ci = R, the ordinary 
Ricci scalar. The invariant of the next order, C2, can be identified with the Gauss-Bonnet (G- 
B) invariant, — 4:R^j_^R^^ + R^^^apR^^""^ ■, which is a topological term as well as a total 
derivative and can be written in 4-d as e^^''^ eap-y&Rfiu "^Rpa '^^ ■ All higher order invariants 
vanish. When we go to 5-d, as is the case we will discuss below, all the £„>3 still vanish as 
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in 4-d but the G-B invariant is no longer a total derivative and its presence will modify the 
results obtained from Einstein gravity in the RS model, altering the equations of motion. 
It is interesting to note that if we demand the absence of ghosts, i.e., terms with no more 
than two derivatives of the metric, then the addition of the G-B piece to the Einstein term is 
the unique modification of the 5-d RS bulk action. What is particularly amazing is that the 
addition of the G-B term to the conventional RS model still allows for a solution with the 
same qualitative structure as is present in the traditional RS model|3- In fact, generalizing 

to D dimensions, Meissner and Olechowski'Sj have shown that RS-like solutions exist with 
the presence of {D — 2)-branes even allowing for all of the non-zero £„ contributions to be 
present in the bulk action. 

Some of the modifications of the RS model due the presence of G-B terms have been 
discussed by other authors {e.g., Refs. |Z1 El El UHl HH 1121 UHl El)- In this paper we are 
interested in how the presence of the G-B term alters the detailed collider phenomenology of 
the RS model. These modifications may occur in many different ways and places. One can 
imagine, e.g., that since the G-B term is of higher order, the triple graviton coupling in the RS 
model^S] may be sensitive to its existence. While this is certainly true, it turns out that this 
is not the most immediate or sensitive way to probe for the existence of G-B contributions to 
the action. As we will see below, the existence of G-B terms will lead to significant changes 
in the masses and SM matter couplings of the Kaluza-Klein gravitons themselves. This 
variation is described by a single new parameter that has a rather restricted range in order 
to avoid tachyons being introduced into the model. In addition, we will demonstrate that 
the anticipated production rate for black holes at colliders can be significantly increased by 
the existence of G-B terms in comparison to the usual RS scenario. For simplicity, we will 
ignore issues having to do with the radion JH] in what follows as this physics may be sensitive 
to the stabilization mechanism[2]. The reader must be careful, however, to insure that no 
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tachyons arise in this sector of the theory when exploring the model parameter space. We 
will also restrict our analysis to the case where the SM fields are confined to the TeV-brane. 
Generalization to the case of bulk SM fields is straightforward [TH] . 

The outline of the paper is as follows: In Section 2 we will present the basic Kaluza- 
Klein formalism and then determine the masses, couplings and wavefunctions for the graviton 
excitations in the RS model in the presence of G-B terms. In Section 3 we will discuss how 
observations of graviton resonances and measurements of their properties at colliders can be 
used to determine the value of the new parameter describing the G-B interactions or place 
stringent upper bound on its value if no deviations from the conventional RS scenario are 
observed. Section 4 contains a discussion of black hole production at the LHC and how it 
will be substantially enhanced by any appreciable G-B terms. We will also show how the 
modifications to the usual RS picture can lead to a significant alteration in the shape and 
parameter dependence of the black hole subprocess cross section. We then summarize and 
conclude. The Appendix outlines the changes to our basic formalism which are necessary if 
graviton brane kinetic terms arising from the usual Einstein action are also present. 

2 Kaluza-Klein Formalism 

The essential ansatz of the 5-d RS scenario is the existence of a slice of warped, Anti-deSitter 
space bounded by two 'branes' which we assume are fixed at the / Z2 orbifold fixed points, 
y = 0, vrrc, termed the Planck and TeV branes, respectivelyPQ. We take the 5-d metric 
describing this setup to be given by the conventional expression 

ds'^ = e-'^^T^^^dx^'dx" -dy\ (1) 

As usual, due to the /Z2 orbifold symmetry we require a = o-{\y\) and, in keeping with the 
RS solution, we expect a = k\y\ with k a dimensionful constant of order the fundamental 
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Planck scale. As first shown in Ref.[7] the inclusion of G-B terms does not alter this basic 
setup. 

The action for the model we consider takes the form 




(2) 



where 



/ 




g —R-Af,+ 



( U"^ /I R _L R tdABCD\ 

— (^/t — ^JXabJ^ + -tlABCD-ti J 



(3) 



describes the bulk with M being the 5-d fundamental Planck scale, A;, the bulk cosmological 
constant and a is a dimensionless constant of unknown sign which we expect to be of order 
unity based on naturalness arguments. Here we note that the upper case Roman indices 
A,B,... run over — 4 while Greek indices will continue to run over — 3. This bulk action 
contains not only the usual Ricci scalar, R, of Einstein gravity but also the Gauss-Bonnet 
quadratic curvature term. Similarly 



describes the two branes with gi being the determinant of the induced metric and Aj the 
associated brane tensions; the Ci describe possible SM fields on the branes. In what follows 
we will assume as usual that the SM fields are all localized on the TeV brane at y = ttVc- For 
simplicity we have not considered potential contributions arising from brane kinetic terms 
for the graviton generated by the Ricci scalars evaluated on the two branes^?]; as we will see 
in the Appendix, such contributions can be included in a relatively straightforward manner. 
G-B like brane terms are automatically absent as the G-B invariant is at most a surface term 
in less than 5-d. 




(4) 
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To go further we insert our metric ansatz into the equations of motion arising from 
the above action; we obtain A^ev = —^pianck as usual together with[T] 

= (5) 

where ' = dy and A{y) is the combination 

A{y) = 6{y) - 6{y - vrrj . (6) 

As in the conventional RS model the solution indeed takes the form a = k\y\ though greater 
than normal care is required to regularize the associated (5-functions. Since 

a" = 2k A , (7) 

using the set of useful relations jHj 

(a')V" = ^k-^'A{y) , (8) 



one then finds 

Apianck = -ATeV = 6kM^[l-^). (9) 

This differs by a factor of 1/3 in the a-dependent term from the work of Kim, Kyae and 
LeejT] due to an error in their application of the proper boundary conditions. In addition, 
using the relation (o"')^" = /c^", we obtain the following explicit expression for k: 



k = k± 



'l±(l ■ ^""^^^'Z' 



2a V V 



1/2 



(10) 
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Note that for a > 0, may in principle take either sign whereas it must be negative if 
a < 0; the conventional RS limit can be obtained by employing the negative root in the 
equation above with A;, negative and taking the limit a ^ 0. 

In order to obtain the equations of motion for the KK excitations we need to go 
beyond this vacuum solution; to this end we consider metric perturbations of the form 



9t,u = e '^"7]^,^ + h 



and employ the traceless-transverse gauge with d^h^y = /i^ = as usual followed by the KK 
decomposition 

V(x,i/) = 5:/i(rj(a;)xn(y). (12) 



such that the Klein-Gordon equation, d^h^^J = —m^h^^J, yields the KK masses. The 
wavefunctions are then seen to satisfyjlUl 



Xn 



d^-4:aa'^d^-8aa'a''dy-4:{a'f{l-4:aa'^)+2a'\l-12aa'^) 



(13) 



where a = a/M^, as obtained by ITU]. Letting Xn = e "^"i^n, the left-hand side of the 
expression above can be simplified using the identities in Eq.(8) to produce 



(14) 



which after some algebra directly leads to the G-B generalization of Eq.(5) in ^HI obtained 
long ago: 



a, 



1-4 



a 

M2 



+ -ne-" l-4^-'' + 4^-'>" = 0- (15) 
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As noted by other authors [3 HUl HH [121 EI, it is straightforward to see that this produces 
the same differential equation in the bulk for the KK eigenf unctions as found in the usual 
RS scenario and thus 

^n(l/) = ^C2(m„e7fc), (16) 

with Cq = Jq^ PnYq being a combination of Bessel functions of order q and where Nn is a 
normalization factor. Following our earlier work^Hl we choose to normalize the KK states 
over the full interval — Trr^ < y < nr^.. For a differential equation with this truncated 
Sturm-Liouville form[TU] 

jO-^Pn = [p{yW + Kw{y)^Pn = , (17) 

the eigenf unctions are to be orthonormalized with respect to the weight function w{y) which 
we can read off in our case from Eq.(15) above: 

My) = e-^'^(l-4a-^ + 8aAA), (18) 

so that 

dy W{y) ^n{y)'4'm{y) = Snm ■ (19) 



Thus, making use of the orbifold symmetry, we obtain the normalization factor 

Nl = 2 1 dye-''y{l - ^a^)^l{y) + Sa—{MO? - e^M^^r^?) ■ (20) 

where e = e~'^^'^''. Neglecting terms suppressed by powers of e, for the massless zero-mode 
graviton this yields the explicit expression 

9 1 / \ 1 / \ / /X 



where we have defined for later purposes the combination 



^ " ■ ^^^^ 

Knowing this normahzation we can directly relate the 5-d parameters M and k to the 
(effective) 4-d Planck mass, Mpi: 

Mi, = ^(l-4«-^)(l + 2Q), (23) 

To obtain the explicit expression for the normalization of the massive KK excitations we 
need to examine the two boundary conditions which are now quite different than in the 
usual RS model. These are obtained by the integration of Eq.(15) above around the two 
branes at y = and Trrc, from this procedure we arrive at 

2 

Til 

V'n(0+)' + O^V(0+) = 0, (24) 

and 

2 

TYl 

M^rt) ' + ne-'^M^rt) = , (25) 

where the eigenfunctions and their derivatives are evaluated on the positive side of both 
branes. These boundary conditions then determine both the wavefunction coefficients /3„ as 
well as the KK masses m„ = Xnke; we obtain 

Cl{Xn) + nXnC2{Xn) = , (26) 

whose roots determine the a;„ and 



Yi{xne) + ^XneY2{xne) ' 
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Note that under most circumstances /?„ ~ and can be safely ignored as in the usual RS 
case; some care, however, needs to be exercised as Q —>■ —1/2. From these expressions we 
can straightforwardly obtain an explicit expression for the normalization factor for the KK 
excitation wavefunctions: 

iVn^o = ^C2(a:„)^(l - Aa^) (1 + 2fi + n'xl) , (28) 

where subleading terms in powers of e have been neglected. The usual RS model is seen to 
correspond to the limit a, f2 — in the expressions above. 

At this point the experienced reader may notice that the modifications to the RS 
expressions above due to the new G-B interactions is strikingly similar in nature to those 
produced by the addition of conventional i?-type kinetic terms on both branes^lTj. This will 
allow one to combine the contributions of both these new physics sources into rather simple 
expressions as we will see in the Appendix below. 




0.0 0.2 0.4 0.6 0.8 1.0 



n 

Figure 1: Imaginary root corresponding to a tachyon in Eq.(26) as a function of positive Q. 
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3 Gravitons at Colliders 



The discussion above tells us that the modifications of the basic RS model due to G-B terms 
can be described by a single parameter, Q. To go further and see this in detail we must 
understand the potential range of this parameter. Since we expect a to be of order unity 
as well as ~ M, one should expect that Q is also of order unity. The arguments given in 
Ref . ^3 suggest that Q cannot be positive as this leads to a tachyon in the KK spectrum. 
To verify this claim we have examined Eq.(26) in detail when i7 > for imaginary roots 
of the form Zn = ±iXn with the results shown in Fig. ^ Here we see that indeed a pair of 
tachyonic roots do exist for i7 > which move off to infinity as we approach the RS limit 
i7 0. We thus agree that solutions with positive values of Q are indeed excluded. A similar 
search for tachyonic roots for i7 < was unsuccessful so we conclude that < to which 
we now restrict ourselves; this implies that the parameter a is also < 0. Eq.(21) provides 
the normalization for the zero-mode graviton wavefunction which must be positive definite 
in order to avoid ghosts. We notice that this requires 4afc^/M^ > — 1 which then implies, 
via the definition Eq.(22) and the negative Q constraint, that Q is confined to the range 

-^<Q<0. (29) 

A short analysis shows that the consideration of the the normalization factor for the gravi- 
ton KK excitations will not improve upon this bound. Thus Q is constrained to be in a 
rather narrow range which significantly increases the predictability of the model. Given the 
definition of Q, Eq.(29) implies that 

A;2 1 

To begin a phenomenological analysis the first question to address is: how does the 
KK graviton spectrum shift as we turn on a non-zero fi? To this end we extract the first few 
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roots, Xn, provided by Eq.(26) and follow their f2-dependence; this is shown in Fig. |21 Here 
we see that the overall Q dependence of the roots is rather weak across the rather restricted 
allowed range for Q. The values of the roots are seen to decrease slightly as Q moves away 
from 0. These small variations in the root values will, however, be of interest to us below 
and will provide an important test of the model. 



15 
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-0.5 -0.4 -0.3 -0.2 -0.1 0.0 

Q 

Figure 2: Behaviour of the first five roots corresponding to the graviton KK masses Xn = 
rrin/ke as functions of the parameter fl. 

Next, we find that the interaction of the KK graviton excitations with the SM fields 
on the TeV brane is given by 

A 

where as usual we define A^^ = Mpie^H]- As in the case of graviton brane kinetic terms, but 
unlike in the usual RS model, the couplings of the KK graviton tower states have become 
level- dependent. In particular, the KK states will always become more weakly coupled as 
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1 + 2fi 
1 + 2VL + n^xi 



1/^ 



(31) 



we go up the KK tower. Furthermore, we see that the couphngs of all the KK states vanish 
as Q —0.5. (It is interesting to note that a very similar rescaling of the conventional RS 
result holds for the case when SM fields are put in bulk.) Fig. IHl shows the rapid decrease in 
the couplings of the first five KK states as Q decreases from 0. The couplings are observed 
to vanish more quickly as one goes farther up the KK tower as expected. 




Q 

Figure 3: Coupling strengths of the first KK graviton states, from top to bottom, in units 
of 1/A^ as functions of the parameter Q. Note that in the RS limit all states have the same 
coupling. 



In addition to specifying Q and the mass of the first KK excitation, mi, to completely 
determine the phenomenology of this model we need to know the ratio c = k/Mpi which 
enters into the expression for the KK widths. In the usual RS model this ratio is bounded by 
the requirement that we don't want potential quantum {i.e., higher curvature) corrections to 
dominate over those from classical Einstein gravity. This requirement is traditionally stated 
as \R\ < M^, which implies that k'^/M'^ < 1/20 for a 5-d Anti-deSitter space; combining with 
Eq.(23) in the a — limit then implies k/Mpi < 1/10 which is the usual result. Under the 
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present circumstances there may be (at least) three objections to this argument. First, since 
we are aheady including G-B terms there is no longer any reason to demand small curvature 
and thus k'^/M'^ may be larger than 1/20. A second objection is that if we do truly want 
small higher-curvature terms for some reason the traditional constraint on fc^/M^ may be too 
weak. To see this we note that for the RS vacuum solution we know that R = —20k'^ (hence, 
the limit above), RabR^^ = 80k^ and RabcdR^^^^ = 40A;^ which implies that the G-B 
invariant has the value 120fc^. These large numerical coefficients may imply that a stronger 
constraint on k'^/M'^ might be necessary. A third observation is that even if k'^/M'^ < 1/20 
holds, for finite a from Eqs.(21) and (23) we find that 



m\i My M3 (1 + 2n) ' 

so that, given the bound on negative Q above, the ratio k/Mpi can become arbitrarily large as 
Q approaches —0.5. While we will restrict ourselves to the conventional RS range of k/Mpi 
in what follows for purposes of comparison, we remind the reader that much larger values of 
this ratio may now be possible. (There are, however, other reasons to believe that this ratio 
may remain as constrained as in the conventional RS model as discussed in Ref.piB ). Given 
a values of and the ratio c = k/ Mpi, we note that the parameter a can now be uniquely 
determined: 

"~ 4c4/3(l + 2^)2/3(1 + ^])l/3 • ^^^) 

With these preliminaries we can now examine how these KK graviton states would 
appear at a collider in comparison to the expectations of the RS model. To this end we 
examine the process e^e~ —>■ ^'^ ^~ in the energy range accessible to the International Linear 
Collider(ILC) and its potential upgrades. A similar analysis can, of course, be performed at 
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Figure 4: Cross section for e+e^ assuming mi = 600 GeV and k/Mpi = 0.05(top) 

or 0.1 (bottom). The usual RS model prediction with = is shown in green while the 
corresponding results in the present model with Q = — 0.2(red) and — 0.4(blue) are also 
displayed. 
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the LHC. To be specific for purposes of demonstration we take mi = 600 GeV and k/Mpi = 
0.05, 0.1; the resuhs of these calculations are displayed in Fig.|3]where we show the cases Q = 
—0.2, —0.4 in comparison to the standard RS model predictions with = 0. Several things 
are immediately apparent: First, as the magnitude of Q increases, the KK states become 
increasingly narrow and their mass splitting is also seen to increase significantly. Second, 
with mi held fixed, the observation of a single graviton resonance and a determination of 
its properties would be insufficient to tell us whether or not G-B terms are present in the 
action. The reason for this is clear: if we just observe a single state we can attribute its width 
solely to the value of A^. We see for example that the width of the first KK is decreased 
as Q decreases from zero by the square of the factor appearing in Fig El The only thing we 
can extract from the first KK width is the value of an effective A^^ which tells us nothing 
about the true value of or Q. To obtain the values of these quantities we need to also 
find the second KK excitation. Fig. E] shows that the mass of the second KK increases 
relative to the first as the magnitude of Q increases; it also gets much more narrow. These 
observations provide our best handles on Q. If we measure the ratio of the first two KK 
masses, -Ri = m2/mi, this will directly tell us the value of Q as shown in Fig. |S| 

We can get a second handle on Q as well as A^^ (and thus k/Mpi through the relation 
"^n = Xn{k/ M pi)At^) by comparing the decay widths of the first two KK states after removing 
'phase space' factors and knowing their respective masses. Since kinematically r„ ~ m^, 
the scaled ratio R2 = r2mf/rim2 would be unity in the usual RS model if we neglect the 
masses of the final state particles. This ratio depends rather strongly on the value of Q 
when G-B terms are present as is shown in Fig. El Recall that as Q approaches —0.5 the 
couplings of the KK states to SM matter on the TeV brane vanishes so that they will no 
longer be produced at colliders. Of course, the ratio of the couphngs of two different KK 
states remains finite. 
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For fixed Q and k/Mpi it is interesting to examine the widths of the graviton KK 
resonances as we go farther up the tower. In the usual RS scenario the widths grow quite 
rapidly r„ ~ and eventually the resonance structure is lost as can be seen in the lower 
panel of Fig.|S| Once G-B terms are present, however, the coupling decreases as m„ increases 
and eventually r„ ~ m„. This means that for appreciable values of Q the KK states are all 
becoming rather narrow as we go further and further up the tower. This can be seen more 
explicitly by examining the ratio r„/m„ as a function of Q for the different KK tower states 
as is shown in Fig. IHl Here we see that for large Q the KK states remain narrow and have 
almost the same value of F/m independent of n. 



B 10"^ 




-0.4 -0.3 -0.2 -0.1 0.0 



n 

Figure 6: The ratio F/m as a function of Q for the first ten KK states. The KK number 
goes up as we go from the bottom to the top of the figure. k/Mpi = 0.05 has been chosen 
for purposes of demonstration. 

One may wonder just how small a value of Q might be measurable. It is clear from 
Figs. lUand El that the ILC will be highly sensitive to small values of Q through the mass 
and width ratios discussed above. Note in particular that in the region near = the Q 
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dependence of R2 is significantly greater than that for Ri. To get a feehng for the level 
of sensitivity we note that a determination of i?i^2 at the level of 1 — 2%, which seems 
rather straightforward at the ILC, will be sensitive to values of Q of order —0.01 to —0.02 or 
less. Thus even small deviations from the standard RS scenario should be observable at the 
ILC. Of course, a detailed analysis accounting for accelerator and detector effects should be 

performed to confirm these conclusions. 

4 Black Hole Production 

The possibility that TeV scale black holes(BH) may be a copious signal for extra dimensions 
at future colliders has been discussed by a number of authors |2II1 E] for both the flat and 
warped background cases. A leading approximation for the subprocess cross-section for the 
production of a BH of mass Mbh is just the geometric BH size [221 IIBI 121] 

a~7ri?2, (34) 

where Rs is the (4 + n)-dimensional Schwarzschild radius corresponding to the mass Mbh- 
The production of BH at the LHC has been studied in detail, including detector effects, in 
Refs. |25| |2n] and there is reason to believe that BH arising from warped and flat backgrounds 
may be experimentally distinguishable [2Z1 due to the RS /Z2 orbifold symmetry or the 
effects of the AdS^ curvature. In our case of interest, we would like to examine how this 
cross section is influenced by including G-B terms in the RS model with n = 1. This can 
be done by following, e.g., the analyses presented in |2H1 1211 1101 IHH 1^21 ESI- ^ simple 
approximate expression for the cross section can be obtained through a modification of the 
fiat space result provided that (Rgke)^ << 1 so that the bulk curvature corrections can be 
neglected; of course we will also need to assume that (e-Rs/27rrc)^ << 1 so that the BH does 
not feel the finite size of the compactified space. (Recall that in the usual fiat space case 
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Rs << 2'n'Rc, where Rc is the compactification radius of the extra dimension.) This is just 
the ordinary Schwarzschild solution in an asymptotically flat 5-d background |28j. We find 
that this approximation works reasonably well for small values of Q but fails when Q nears 
its lower limit Q = —0.5 since the AdS curvature is larger in that case. To obtain the full 
expression |30j these bulk curvature terms can no longer be neglected and we must consider 
the BH to be embedded in (an infinite!) AdS^; note that this expression still assumes that 
Z = (e-Rs/2vrrc)^ << 1, which remains as an approximation. We will see that this inequality 
is reasonably well satisfied below. Once we adjust for the definitions of the fundamental 
parameters as given in the action above, one obtains 



TC 

a 



2/3M2 



1/2- 



1 + 1 - 4/3(2a - 7) 



(35) 



where 



Mbh 

^ = 3^' ^''^ 

and where is the 'warped-down' fundamental scale, M^, = Me. This expression is not 
very transparent so let us expand the square root to first order; we then find 



. Mbh 
a ' 



37rM3 



1 — Gvr^a 



BH 



(37) 



which is the leading order term in the bulk curvature expansion; the usual RS result in the 
flat space approximation is obtained by setting a = in this expression| 



To evaluate the cross section we note that Eq.(33) tells us that a is determined 
provided we know both Q and c = k/Mpi. Note that for k/Mpi flxed, a becomes quite 
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large (and negative) as f2 ^ —0.5. Given the large prefactor of Gvr^ in the approximate 
expression above we might expect that G-B terms can be numerically quite important even 
for small values of Q. It is interesting to note that this cross section is bounded from below 
as can be seen from the exact expression, i.e., a ^ —l-aajMl, independent of the BH mass. 
(Remember that a is negative here.) In obtaining our numerical results we will make use of 
the exact expression above. 

Fig. [7| shows the influence of the G-B terms on the BH production cross-section at the 
LHC; for definiteness we have assumed M^, = 1.5 TeV and kjMpi = 0.1 but there is nothing 
special about these particular values. We see, as expected, that a non-zero fl can lead to 
a significant increase in the BH production cross section; for Q = —0.1 (—0.4) this increase 
is seen to be approximately one(two) order(s) of magnitude in comparison to the usual RS 
scenario. Here we may need to worry that the approximation that the effects of the finite 
radius of compactification can no longer be neglected. The inequality Z = {eRs/2'n-rc)^ << 1 
is not so obviously satisfied here since the value of Rs is growing so large. To test this let A 
be the ratio of the subprocess cross-section in the presence of G-B terms to that in the RS 
model. Then the quantity Z can be expressed as 



so that, with k'^/M^ = 1/20 and Mbh = 5M^, for purposes of demonstration, we arrive at 
Z ~ 1.7(A/100) X 10~^ which is significant only for very large A. Thus we see that for almost 
all of our parameter space Z << 1 remains valid. 

Although we have seen a substantial increase in the BH cross section relative to the 
RS model due to G-B terms, the overall normalization of the cross section remains uncertain 
due to the complexities of BH formation by terms of order unity as has been discussed by 




(38) 
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Figure 7: (Top) Black hole production cross section at the LHC assuming = 1.5 
TeV and k/Mpi = 0.1 as described in the text. The lowest dotted curve corresponds 
to the RS model prediction with f2 = 0; the subsequently higher curves corresponds to 
fl = —0.02, —0.1, —0.2, —0.3 and —0.4, respectively in the G-B case. (Bottom) Same as in 
the top panel but now for a fixed minimum black hole mass of 6 TeV as a function of Q. 
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many authors [221 123 1^- While overall factors may be uncertain we might expect the various 
scaling laws to remain at least approximately valid. In the case of flat extra dimensions, 
or in the RS model in the {RskeY << 1 limit, the value oi n = d — A uniquely determines 
the simple dependence of the cross section on the ratio Mbh/M^: aM^ ~ [Mbh / M^]'^^^"~^^\ 
Here, for the case of a single warped extra dimension, the G-B terms conspire to modify the 
shape of the cross section as a function of Mbh/M^, in a rather unique manner. This can 
even be seen by examining the approximate expression Eq.(37): the subprocess cross section 
is generally no longer proportional to a single power of Mbh/M*- In fact, as the G-B terms 
come to dominate, the subprocess cross section becomes essentially independent of Mbh as 
shown in Fig. |S1 where the important feature to notice is the shape of the distribution and 
not the overall magnitude. Note that this effect happens quite rapidly as a non-zero value 
of Vt is turned on. This rather distinctive behavior may be observable in BH production at 
the LHC once the parton distribution luminosity is unfolded and will help to pin down the 
existence of G-B interactions. It should then be possible to correlate shifts in the production 
properties of BH with those of graviton resonances to precisely determine the G-B parameter 
VL. 

5 Conclusions 

Higher-dimensional curvature corrections to Einstein gravity are to be expected on rather 
general grounds; the observation of any effect arising from such terms will tell us important 
information about the ultraviolet completion of this theory. The requirements of unitarity 
and the absence of ghosts in the 5-d Randall-Sundrum model uniquely determine the form 
such terms may take and eliminates any derivatives higher than second in the equations 
of motion. Consequentially, the Gauss-Bonnet invariant is singled out as the only possible 
new non-trivial addition to the bulk action for gravity beyond the usual Ricci scalar. The 
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Figure 8: The BH sub-process cross section as a function of the ratio Mbh/^* for different 
values of fi. From bottom to top the curves correspond to f2 = 0, —0.03, —0.1, —0.2, —0.3 
and —0.4, respectively. 

existence of G-B terms in the action has been previously shown to leave the qualitative 
features of the RS model invariant but the quantitative details are expected to be altered, 
perhaps significantly. 

In this paper we have begun to examine the phenomenological implications of the 
existence of additional G-B terms in the RS bulk action. Two of the dominant signatures 
for the RS model are the production of TeV-scale graviton resonances and the production 
of black holes at future colliders. This paper has shown that both of these processes can be 
significantly altered if G-B terms are present in any appreciable amount. If no deviations 
from the classic RS model predictions are observed, we have shown that stringent bounds 
can be placed on the parameter Vt which describes the relative strength of the G-B terms. 

The first step in this analysis was to determine how the graviton KK masses, couplings 
and wavefunctions are modified by the presence of the G-B terms. Once this was done we 
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analyzed the production properties of the KK graviton resonances at the ILC; we showed 
the mass spectrum and decay widths of these resonances were substantially different than in 
the RS case over most of the G-B model parameter space. It was demonstrated that the KK 
mass and width ratio measurements at the ILC can be highly sensitive to the presence of 
G-B terms and can be used to precisely determine their strength or to place strong bounds 
on their existence if no deviations from the RS model are observed. 

Next we examined the production of BH at the LHC in the RS model with G-B terms. 
We found that not only is the production cross section significantly enhanced in comparison 
to the conventional RS model but its parametric dependence on the ratio Mbh/M^ can be 
drastically altered even for small values of Q. Any deviations observed in the production of 
graviton resonances at the LHC or ILC can then be directly correlated with modifications 
to BH production to determine the value of Q. 

Hopefully signals of warped extra dimensions will be observed at future colliders and 
we may be ultimately be able to probe the ultraviolet completion of quantum gravity. 
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Appendix 

In this Appendix we consider a modification of the action given by Eqs.(3) and (4) above 
to include graviton brane kinetic terms generated by the 4-d Ricci curvature scalars on the 
TeV and Planck branes. Such terms have been suggested in the case of G-B extensions of 
the RS model^l] to help assist removal of the tachyon field present |T3"^ when a > 0. Using 
the notation above as well as in Ref . [T7j , this amounts to adding to the action a term of the 
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form 

M3 



Sx4^ — yg^Tc^f {y) + gnrcS {y - vrrj J R^'^^ , (39) 

where R^^^ is the induced 4-d curvature and (?o,7r are dimensionless constants that we expect 
to be of order unity. Since these are boundary terms they do not modify the bulk equation 
for the KK wavefunctions which remain unaltered from the traditional RS case as discussed 
above. Defining the combinations 

70,^ = 9o,7T^ , (40) 

as in our earlier work|17j. we find that the only change induced by these new brane terms is 
to modify the definition of w{y) given in Eq.(18) above to 

w{y) = e-'^ (l - + 8a^A + ^5(y) + ^5(y - vrrj j , (41) 

with the subsequent appropriate changes in the boundary conditions: 

2 

TTI 

^n(0+)' + fio^^(0+) = 0, (42) 



and 



2 

Ttl 

Mnrt) ' + r].e- V^„(7rr+) = , (43) 
where now we now define the combinations 

The generalizations of the other equations in the main text can now be read off directly, e.g., 

Cl(x„) + Vt^XnC2{Xn) = , (45) 
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and 

Ji{xne) + QoXneJ2{ 

Fi(x„e) + VtQXneY2{xne) ' 

determine the masses and wavef unctions of the KK states. A short analysis shows that there 
will exist parameter ranges that are tachyon free. The collider and black hole analyses will 
now require a detailed exploration of this VLq — VLt^ parameter space which is beyond the scope 
of the present paper. 
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